HLRZ 65/93 



INSTANTONS AND SURFACE TENSION 
AT A FIRST-ORDER TRANSITION 



Sourendu Gupta 
HLRZ, c/o KFA Jiilich, D-5170 Jiilich, Germany. 



ABSTRACT 

We study the dynamics of the first order phase transition in the two 
dimensional 15-state Potts model, both at and off equilibrium. We find 
that phase changes take place through nucleation in both cases, and finite 
volume effects are described well through an instanton computation. Thus 
a dynamical measurement of the surface tension is possible. We find that 
the order-disorder surface tension is compatible with perfect wetting. An 
accurate treatment of fluctuations about the instanton solution is seen to be 
of great importance. 



HLRZ 65/93 
February 2008 



First order phase transitions, i.e., phase coexistence points, have recently been 
subjected to intense analysis. Such systems are characterised by many different 
dimensionful quantities. It is useful to divide these into two classes of observables. 
The first pertains to properties of the pure phases. Such properties are obtained, 
as usual, by taking derivatives of the (extensive) free energy. These derivatives, or 
cumulants, are extracted through finite-size scaling. In recent years a full theory of 
such scaling has been developed [1,2] and tested [1-3]. The second class concerns 
coexistence. Foremost among such variables is the surface tension — the leading non- 
extensive part of the free energy. Most theories of the (canonical) dynamics at the 
phase transition involve the surface tension. The measurement of this quantity is 
usually approached through detailed investigations of the static system at the phase 
transition [4]. 

However, much of the recent attention enjoyed by first-order phase transitions 
is due to the interesting dynamics of the transition. This is expected to be due 
to nucleation. We perform a careful analysis of finite size effects in the equilibrium 
dynamics and compare our observations with an instanton-based computation [5] . We 
also test nucleation theory through a non-equilibrium process — hysteresis. Certain 
scaling laws for this have been proposed recently [6] ; we verify them for the first time. 
These two tests provide the justification for the use of dynamical techniques for the 
measurement of the surface tension. 

Our numerical work is done with a simple model — the two dimensional 15 state 
Potts model. This model is defined through the partition function and Hamiltonian 

Z = J2^M-PH), H = J] 1-5,^,,, (1) 

Wr} {rr') 

where the spin sitting at the site r of a (square) lattice can take one of 15 values 
and the angular brackets denote nearest neighbours. At j3t = log(vT5 -|- 1), ordered 
and disordered phases coexist in the thermodynamic limit. The two phases can be 
distinguished either through a singlet magnetisation or the internal energy density 

e = ^{H). (2) 

where the lattice size is L?. In the thermodynamic limit, the internal energies in the 
ordered and disorder phases are, respectively, cq = 1.779 and = 0.737. 

Recent analysis of correlation lengths in Potts models [7] have yielded an exact 
result for a spin-spin correlation length. This has been identified with the correlation 
length in the disordered phase ^d- A duality argument has been used [8] to relate the 
order-order surface tension, (Too-, to this correlation length — 

Goo = (3) 
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The perfect wetting conjecture would then imply the relation 



2a-od = (Too = 



(4) 



between the order-disorder surface tension, cTod, and the rest of these quantities. 
Although proven only in the limit — > oo, it is strongly suspected that perfect 
wetting holds for two dimensional Potts models for all q> A. From the formulae of 
[7] one finds = 4.18 for q = 15; implying 2aod = 0.239. Thus a quantitative test 
of nucleation theory should yield this value for aod- Alternately, the argument can 
be turned around and the measurement of the single quantity, aod, can be used to 
check the perfect wetting conjecture. 

One of the dynamical methods we use is a finite-size scaling of the exponential 
autocorrelation time, r, of the energy density. Through a study of the 10-state 
Potts model, numerical evidence was presented in [9] that this autocorrelation time 
is determined by the tunnelling rate between the coexisting ordered and disordered 
phases. As a result, the autocorrelation time can be identified with the outcome of 
an instanton computation [5] giving 



In our case, of course, d = 2. The exponential factor is the saddle-point result; and 
the L-dependence of the prefactor is obtained from a one-loop computation of the 
determinant of the fluctuations around the saddle point. The test of the instanton 
computation is in the L-dependence of the measured values of r. Note that Eq. (5) 
describes dynamics in equilibrium. 

We also study a particular example of off-equilibrium dynamics, that of hystere- 
sis. The coupling f3 is cyclically varied about the critical coupling Pt with a frequency 
u) and an amplitude A/3. Hysteresis occurs as e switches between the values Co and 
Cd- The area of the hysteresis loop in the energy density, A{lj) is studied as a func- 
tion of u>. The system has a real-time excitation with a 'mass' given by Eq. (5). 
The measurement of A{(jj) is really a 'line-shape' measurement. When A/3 — > 0, the 
peak should be at t~^{L) and the shape should be approximately Lorentzian. This 
cannot be converted to an useful test because the functional form of the finite A/3 
corrections is not known. 

We define the coercive coupling SPc through the fact that hysteresis loops reach 
the value (cq -|- e^) /2 when the coupling is Pt i SPc- Here the nucleation rate becomes 
larger than 1/a;, and the probability of a flip into the stable phase exceeds 1/2. This 
argument was presented in [6] and was developed into the scaling law 



T-\L) = aL-'^/'^exp{-2aodL'^~^)- 



(5) 



A{u) - (Ina;)-^/^'^-^) where a; ^ 0. 



(6) 
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The coercive coupling also obeys the same scaling law. Verification of this relation 
thus furnishes a test of nucleation theory. Furthermore, since the tunnelling rate is 
finite for any finite lattice, A{u>) is zero at some non-zero frequency Q, and we have 
the relation 

J)(L) = T-^(L) where ^(11) = 0. (7) 

Thus the scaling of Q with L is again given by Eq. (5), and constitutes yet another 
test of nucleation theory. Note that Eq. (6) refers to a slow non-equilibrium situation. 

TABLE 1. 

Run parameters for the two-dimensional 15-state Potts model. We show the 
lattice sizes L, pseudo-critical couplings /3t{L), statistics used for the determination 
of autocorrelation times, Nmn, and the hysteresis parameters Pm, and Sp. 



L 


PAL) 


■^run 


set 1 


set 2 


set 3 


Pm 


A/3 


5P 


Pm 


A/3 


5P 


Pm 


Ap 


6p 


8 


1.546 


1 X lO'^ 


1.65 


0.200 


0.010 














12 


1.5661 


1 X 10*^ 


1.65 


0.100 


0.005 


1.61 


0.060 


0.005 


1.6470 


0.1620 


0.0135 


16 


1.5723 


2 X 10^ 


1.61 


0.080 


0.005 


1.61 


0.060 


0.005 


1.6260 


0.108 


0.009 


20 


1.5772 


4 X 10^ 


1.61 


0.060 


0.005 














24 


1.5798 


5 X 10^ 




















30 


1.5811 


4 X 10^ 





















For Potts models at phase coexistence it was observed [9] that both local and 
Swendson-Wang dynamics are dominated by tunnellings, and that the exponential 
autcorrelation times, with changing P and L, are related by a constant. In view of 
this, all our simulations were performed with the latter algorithm. Autocorrelation 
times were measured at the pseudo-critical couplings, /3t(L), defined by the maximum 
of the specific heat. For L < 20 the values of Pt{L) and r were obtained in [10]; 
this work verifies these measurements. Hysteresis was induced by cyclically changing 
P from a maximum of Pm down by an amount A/3 and back, in discrete steps of 
5/3, running N sweeps of a Swendson-Wang update at each coupling. The runs were 
started by first thermalising a system at Pm with 2 x 10^ cluster updates. Then for 
each (and fixed values of the other parameters) we ran through 200 hysteresis 
cycles. Since fairly large values of N had to be used, this was by far the most CPU- 
intensive part of these computations. The run parameters are shown in Table 1. 
Three difi'erent sets of parameters were used in order to check that the scaling law of 
Eq. (6) and the extrapolated values of il[L) were independent of A/3. 



3 



For each hysteresis loop, we defined the area by the sum 

A = 5/3^(-l)^e-, (8) 

i 

where i labels each of the values of f3 in the cycle, the bar above denotes averaging 
over the sweeps performed at that f3. The value of D was set to be equal to in 
that half of the cycle with decreasing /3, and 1 in the other. The averages and errors 
were obtained by jack-knife estimators over all the cycles. 

A second measurement was of the cyclic response function 

ct = {{e^ - {em, (9) 

where the angular brackets denote averages over all measurements performed at a 
given coupling in the hysteresis cycle. This response function peaks twice during a 
cycle, at ±5/3c, and allows us to extract SPc by searching for the maximum of c^. 
This procedure is the dynamical analogue of defining the transition coupling on a 
finite lattice, Pt{L) by the peak in usual response function Cy 

The average and error were again estimated by a jack-knife procedure. These 
measurements of SPc have larger relative errors than A. This is due to two reasons. 
The first is intrinsic. Since tunnellings occur at random times, there is a cycle to cycle 
variation in the coupling at which tunnellings occur. The other error is related to the 
statistics. The identification of the coercive coupling depends on the relative heights 
of the two peaks, and is subject to fairly strong errors. Due to these uncertainties, 
we decided not to use this quantity for our scaling tests. 

The frequency oj should be identified with 27r5p/NAp. However, for each L, 
since SP as well as A/? are fixed inside each set, these factors are not important when 
trying to check the scaling with u. Furthermore, when comparing different values of 
L, we are interested in time scales expressed directly in sweeps. Hence we shall use 
the convention cv = 1/N. This is only a matter of convenience. When necessary, one 
should use the full definition of w. 

For each L and a set of uj at fixed A/9 and 6(3, we tried to fit the data on 
A{lj, A/3) to a Lorentzian. We found that this description improves as A(3 decreases. 
This is illustrated in Fig. 1 for the L = 16 lattice. At low frequencies, where the data 
deviates from the Lorentzian shape, A{u) is independent of Aj3. For large N, when 
the data deviate from the Lorentzian shape, we fit a form 

A{uj) = a-b/\ogN. (10) 

We found extremely good fits to this form and could certainly rule out any power-law 
behaviour. The data and fits are shown in Fig. 2. The fitted parameters a and b give 
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1/N 

FIGURE 1 

Hysteresis loop areas A{(jj) as functions of 1/N for the 15-state Potts model on 
16^ lattices. The data are for sets 1 (squares), 2 (circles) and 3 (triangles) of Table 
1. The lines show the best Lorentzian fits. 

estimates of 0, — the frequency at which the loop area vanishes. Errors on were 
estimated from the covariance matrix between these parameters. We found good 
agreement between these values obtained indirectly and the direct measurements of 
t(L). The scaling of Q,~^{L) with L is consistent with Eq. (5), but does not provide 
a very stringent test. The direct measurements of t(L) are, of course, more accurate. 

The direct measurements of t(L) were obtained by constructing the autocorre- 
lation function and fitting its long-distance form to an exponential. As a cross check, 
we measured local masses and looked for plateaus as a signal that a single mass fit 
over a given range was reasonable. The measurement procedure remains the same 
as in [9]. In all cases the fits were performed over a range which turned out to be 
T < t < 4:T. The errors on t(L) were, of course, reflections of the errors on the 
autocorrelation function. These were obtained as the dispersion between jack-knife 
blocks. We varied the number of jack-knife blocks between 5 and 25. We took the 
lack of sensitivity of the means and errors to the number of blocks as an indication 
that our error estimates are reliable. 

In order to test Eq. (5) and measure aod, we fitted the data on t(L) to the form 

log(T/L) = cL + c'. (11) 
Note that aod is given by c/2. On the left hand side of Eq. (11), the division by 



^.04 
3^ 




0. 18 
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FIGURE 2 

Hysteresis loop areas A{(jj) as functions of l/logA?" for the 15-state Potts model 
for L = 8 (filled circles), 12 (filled circles), 16 (filled squares) and 20 (open squares). 
The lines show the best fits. 

L takes care of the effects of fluctuations around the instanton solution. It turns 
out that this term in the fit is quite crucial. An attempt to perform the fit without 
this factor was completely unsuccessful; values obtained increased by almost an 
order of magnitude. In principle, one could perform a three-parameter fit, leaving the 
power of L in the pre-exponential factor to be determined by the data. Unfortunately 
this requires more lattice sizes than we had in this study. Our fits gave 

2(7od = 0.263 ± 0.009 (x^ = 1.0/2 dof, 16 < L < 30). (12) 

This result is compatible with perfect wetting. 

Although the difference between the perfect wetting result for Uod and our mea- 
surement is not statistically significant at the 3a level, we believe it deserves comment. 
We find it difficult to regard seriously the possibility that perfect wetting begins to 
break down when q drops to a number close to 15. More likely is that Eq. (5) has 
to be supplemented with a higher loop computation. It has been argued [5] that a 
loop- wise expansion of the pre-exponential factor yields a power series in L~^/'^. It 
is a reasonable guess that the two-loop term is marginally important for the lattice 
sizes we have worked with. Then our observations would imply that the coefficient 
of the term in is positive. A computation of this term would certainly be useful. 
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FIGURE 3 

We show the scaHng of 1/0 (open circles) and directly measured values of the 
autocorrelation time at Pt{L), r, (filled circles) against the lattice size L. The lines 
show the best fits of the form shown in Eq. (13). The values of have been 
multiplied by 2 for visibility. 

Finally we comment on previous numerical tests of Eqs. (5) and (6). A high- 
statistics study of the 10-state Potts model [9] had established that the autocorre- 
lation time in equilibrium was determined by the tunnelling phenomenon. However, 
this study had not been able to observe even the dominant exponential behaviour in 
Eq. (5). It was conjectured there that L/^ values used there were too small. The 
recent work of [7,11] shows that this is indeed correct. In that study the largest 
values of L used were about 3^, whereas this study uses between 4^ and 7^. 

Earlier studies of the scaling of A{lij) with uj had parametrised the variation 
by power laws. This is presumably correct for some systems, but the arguments 
of [6] must hold whenever the dominant dynamical mechanism is nucleation and 
tunnelling. Magnetic hysteresis in the Ising model, or the one-component (f)'^ theory 
should therefore be described by Eq. (6). The contradictory results of [12] were 
obtained with values of u much smaller than the ones we use. The implication is 
that Eq. (6) is not applicable to these. 

We summarise the main results obtained in this study. The scaling law of Eq. (6) 
for the frequency dependence of hysteresis loop areas is found to hold extremely well 
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over three decades in frequency, and for a variety of lattice sizes. This is strong evi- 
dence that in the particular non-equilibrium situation at a first-order phase transition 
exemplified by hysteresis, the dynamics is of nucleation. Furthermore, the expression 
in Eq. (5) is found to describe the finite-size scaling of the autocorrelation times, 
showing that an instanton based description of the equilibrium dynamics is valid. A 
proper treatment of fluctuations around the instanton is observed to be crucial for 
the description of the data. The dynamics then allows the extraction of the surface 
tension. For the 15-state Potts model we find that perfect wetting holds. A statis- 
tically insignificant discrepancy can be attributed to the neglect of two-loop terms 
in the treatment of the fiuctuation determinant. It should be emphasised that this 
makes the present computation one of the most accurate measurements of a surface 
tension to date. 



8 



REFERENCES 

1) C. Borgs and R. Kotecky, J. Stat. Phys. 61 (1990) 79; 

C. Borgs, R. Kotecky and S. Miracle-Sole, J. Stat. Phys. 62 (1991) 529; 
C. Borgs and W. Janke, Phys. Rev. Lett. 68 (1992) 1738. 

2) S. Gupta, A. Irback and M. Olsson, preprint HLRZ 23/93 and LU-TP-93-6, 
NucL Phys. B, in press. 

3) A. Billoire, R. Lacaze and A. Morel, Nud. Phys. B370 (1992) 773; 
B. Berg, A. Billoire and T. Neuhaus, Saclay preprint SPhT-92/120. 

4) K. Binder, Z. Phys. B 43 (1981) 119; 

K. Kajantie, L. Karkkainen and K. Rummukainen, Phys. Lett. B223 (1989) 213; 
J. Potvin and C. Rebbi, Phys. Rev. Lett. 62 (1989) 3062; 
K. Jansen et al, Nucl. Phys. B322 (1989) 693. 

5) J. C. Niel and J. Zinn-Justin, Nucl. Phys. B280 [FS18] (1987) 355. 

6) D. Dhar and P. B. Thomas, preprint TIFR/TH/92-32. 

7) E. Buffenoir and S. Wallon, Saclay preprint, SPhT/92-077. 

8) C. Borgs and W. Janke, preprint HLRZ 54/92 and FUB-HEP 13/92. 

9) A. Billoire, R. Lacaze, A. Morel, S. Gupta, A. Irback and B. Petersson, Nucl. 
Phys. B358 (1991) 231. 

10) M. Ohlsson, B. Sc. Thesis, University of Lund, unpublished. 

11) S. Gupta and A. Irback, Phys. Lett. B286 (1992) 112. 

12) M. Rao, H. R. Krishnamurthy and R. Pandit, Phys. Rev. B42 (1990) 856; 
W. S. Lo and R. A. Pelcovits, Phys. Rev. A42 (1990) 7471. 



9 



